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7 The detailed analysis of the structural variations of three GeC>2 and SiC>2 polymorphs 
s (a-quartz, a-cristobalite, and rutile) under compression and expansion pressure is re- 

9 ported. First-principles total-energy calculations reveal that the rutile structure is the 

10 most stable phase among the phases of GeC>2, while SiO"2 preferentially forms quartz. 

11 GeO"4 tetrahedras of quartz and cristobalite GeO"2 phases at the equilibrium volume are 

12 more significantly distorted than those of Si02- Moreover, in the case of quartz Ge02 

13 and cristobalite Ge02, all O-Ge-0 bond angles vary when the volume of the Ge02 bulk 



changes from the equilibrium point, which causes further deformation of tetrahedra. In 



15 contrast, the tilt angle formed by Si-O-Si in Si02 markedly changes. This flexibility of 

OO ' 

Q\ ! i6 the O-Ge-0 bonds reduces the stress at the Ge/Ge0 2 interface due to the lattice-constant 
G\ ■ 

17 mismatch and results in the low defective interface observed in the experiments [Matsub- 
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is ara et al: Appl. Phys. Lett. 93 (2008) 032104; Hosoi et al: Appl. Phys. Lett. 94 
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21 1. Introduction 

22 Ge has recently attracted increasing attention for future advanced complementary metal 

23 oxide semiconductor (CMOS) device structures owing to its high intrinsic carrier mobil- 

24 ity as it is becoming increasingly difficult to enhance the performance of CMOS devices 

25 through scaling based on conventional Si-based techniques. The key issue to be resolved in 

26 advanced Ge-based devices is the formation of gate stacks with superior interface proper- 

27 ties. Research activity related to Ge0 2 is steadily increasing, and it has been reported that 

28 Ge is well passivated by Ge0 2 by conventional dry oxidation without any hydrogen pas- 

29 sivation treatment. 1,2 On the theoretical side, Houssa et al. claimed that the viscoelastic 

30 properties of Ge02 lead to a low interface defect density at the Ge/Ge02 interface af- 

31 ter performing a calculation using a modified Maxwell model. 3 In our previous study, 

32 we examined the oxidation mechanism of crystalline Ge0 2 and Ge(100) interfaces by a 

33 first-principles total energy calculation following the Si(100) oxidation process proposed 

34 by Kageshima and Shiraishi, 4 and we found that Ge atom emission, which deteriorates 

35 the Ge/Ge0 2 interface, hardly occurs during the oxidation process of Ge(100). 5 To inves- 

36 tigate the low probability of Ge atom emission, the mechanism accounting for the release 

37 of interface stress should be clarified. 

38 From the viewpoint of phase transition, Si0 2 occurs in many different forms. At ambient 

39 temperature and pressure, the ground-state structure for Si0 2 is a-quartz (q-Si0 2 ). Si0 2 

40 forms a rutile structure (r-Si0 2 ) under a pressure above 2 GPa and it is transformed 

41 into a-cristobalite (c-Si0 2 ) at a high temperature. On the other hand, there are two 

42 stable polymorphs of Ge0 2 at normal pressures: the low temperature form has the rutile 

43 structure (r-Ge0 2 ), and Ge0 2 undergoes a smooth transformation to a-quartz (q-Ge0 2 ) 

44 at T 1300 K. An a-cristobalite structure (c-Ge0 2 ) has been identified after the long- 

45 time heating of Ge0 2 glass at 873 K. 6 The investigation of the c-Ge0 2 phase is important 

46 because Ge substrates are typically subjected to dry oxidation in an 2 ambient at 623-823 

47 K to form a Ge0 2 layer in recent experiments. 2 ' 7 

48 Here, we investigate the structural properties of q-, c-, and r-Ge0 2 by first-principles 

49 total energy calculations. The structural properties of Si0 2 are also examined for compar- 
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50 ison. It was found that the rutile structure of Ge0 2 is the most stable structure, whereas 

51 Si0 2 preferentially forms the quartz structure. We then investigate the variations of the 

52 O-Ge-0 (O-Si-O) bond angles of quartz and cristobalite phases with respect to the vol- 

53 ume since local pressure is induced at the semiconductor /oxide interface. Although the 

54 variations of the atomic structures of q-GeC>2, r-GeC>2, q-Si02, c-Si02, and r-Si02 under 

55 pressure have been examined by both experimental and theoretical studies, 8-20 no study 

56 reports the relationship between their bond angles and pressure of these six oxides in 

57 the same treatment of computational code or experimental facility. We note that the 

58 lattice constant (bulk modulus) of q-Si02 reported by the first-principles calculation 15 is 

59 smaller (larger) by more than 3 % (34 %) than the one reported by another first-principles 
eo study. 16 Moreover, it is reported that a dioxide forms a cristobalite structure before the 
ei atom emission as well as at the initial stage of oxidation. 4 The pressure-dependent behav- 

62 ior of c-Ge0 2 , which has never been explored to the best of our knowledge, is a subject of 

63 intense research to clarify the relaxation mechanism of the interface stress. Therefore, it 

64 is of importance that the uniform theoretical treatment facilitates systematic comparisons 

65 and the identification of trends among these six oxides. Our finding is that the O-Ge-0 
ee bond angles change significantly under pressure, while the tilt angle and Si-0 bond length 
67 vary in the case of Si0 2 . The variation of the bonding network of Ge0 2 exhibits com- 
es pletely different characteristics from that of Si0 2 . The metallic properties of Ge provide 

69 a qualitative understanding of not only the difference between the ground-state phases of 

70 Ge0 2 and Si0 2 but also the variations of the bond angles under pressure. 

71 The rest of this paper is organized as follows. In § 2, we describe the computational 

72 techniques used in this study. In § 3, we present the main results and a discussion of the 

73 structural parameters and properties. Finally, a brief summary is given in § 4. 

74 2. Computational Techniques 

75 The structures of q-, c-, and r-Ge0 2 are hexagonal, tetragonal, and tetragonal with 

76 three, four, and two Ge02 molecules per unit cell, respectively. Ge atoms in q- and c- 

77 Ge0 2 are surrounded by four oxygen atoms, while each Ge atom in r-Ge0 2 is surrounded 

78 by six oxygen atoms with distorted octahedral coordination as shown in Fig. 1. The 
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79 calculations are performed within the local density approximation 21 of density functional 

so theory 22,23 using the real-space finite-difference approach 24-28 and the norm-conserving 

si pseudopotentials 29 of Troullier and Martins 30 in the Kleinman-Bylander representation. 31 

82 The grid spacing was set at 0.25 bohr, and a denser grid spacing of 0.083 bohr in the 

83 vicinity of nuclei with the augmentation of double-grid points 27 for each Ge0 2 polymorph. 

84 We took 4x4x4, 4x4x3, and 4x4x6 Appoint grids in the Brillouin zone for q-, c-, 

85 and r-Ge0 2 , respectively. The optimal lattice parameters and internal atomic coordinates 
se were determined by minimization of the total energy using calculated forces, with a force 
a? tolerance of F max < 1.0 mH/bohr. The same computational procedures were applied for 

88 Si02 polymorphs. 

89 3. Results and Discussion 

90 Figures 2 (a) and (b) show the total energy per molecular unit as a function of volume 

91 for the three structures of Ge0 2 and Si0 2 , respectively. The zeros of the energy scales 

92 are rutile for Ge0 2 and quartz for Si02- It was found that the zero-temperature phase in 

93 Ge0 2 has the rutile structure while that in Si0 2 exhibits the quartz structure, which is in 

94 good agreement with a previous report. 18,32 Sn, which is also a group IV element similarly 

95 to Si and Ge, is a metal and its oxide crystallizes in the rutile structure under ambient 

96 conditions. 33 Ge0 2 forms the sixfold-coordinated rutile structure more preferentially than 

97 Si0 2 because Ge is between Si and Sn in the periodic table. Tables I and II show the 

98 calculated lattice constants, bond lengths l i: bond angles and tilt angles 5 of Ge0 2 and 

99 Si0 2 . The other calculated and experimental results are also given. 9, 10,20,34-37 There are 

100 two distinct Ge-0 (Si-O) bond lengths in Ge0 4 (Si0 4 ) tetrahedra. In addition, 0j exactly 

101 corresponds to the O-Ge-0 (O-Si-O) bond angle, and the tilt angle is related to the Ge-O- 

102 Ge (Si-O-Si) bond angle. 38,39 The agreement between our results and experimental results 

103 for Si0 2 is excellent for the structural parameters. The lattice constants of q- and c-Ge0 2 

104 are slightly underestimated in both the theoretical calculations: this underestimation is 

105 caused by the use of the local density approximation, 21 and the parameters obtained by 
we the theoretical calculations agree well. The deviations of the O-Ge-0 bond angles from 
107 the ideal tetrahedral angle (109.5°) are larger than those for O-Si-O, resulting in distorted 
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108 Ge0 4 tetrahedra. 

109 Thermally oxidized Ge0 2 on a Ge substrate has been found to mainly form the fourfold- 
no coordinated structure by X-ray photoelectron spectroscopy, 40 and crystalline q- and c- 
111 Si02 are formed on a Si(100) surface according to the oxidation model of the Si(100) 
n2 surface. 4 The compressive in-plane stress at the Ge/Ge02 (Si/Si0 2 ) interface is induced 
n3 by the lattice mismatch between Ge (Si) and its oxide. Therefore, we particularly focus on 
ii4 the bond structures of crystalline q- and c-Ge0 2 assuming that the oxidation mechanism 
us of Ge is the same as that of Si. Figures 3, 4, and 5 show the variations of the Ge-0 
lie (Si-O) bond lengths, the O-Ge-0 (O-Si-O) bond angles, and the tilt angles with respect 
ii7 to the volume, respectively. In the case of the rutile phases, we depict the variations of 
us the Ge-O-Ge (Si-O-Si) bond angles instead of the tilt angles since the rutile structure has 
n9 a higher symmetry than the others. The models at elevated pressures do not exhibit the 

120 amorphous phase or transform into another phase because the calculated O-Ge-0 and 

121 O-Si-0 bond angles show no indication that the tetrahedra become significantly more 

122 irregular or distorted. Although the lattice constants are varied by increments of 1%, 

123 the variations of the bond lengths are less than ~0.1%. This indicates that the bond 

124 angles play a predominant role in the compression or expansion of the fourfold oxides. 

125 Note that the bond lengths in the sixfold oxides change under pressure due to the higher 

126 symmetry of the O-Ge-0 (O-Si-O) bond angles as well as the Ge-O-Ge (Si-O-Si) ones. 

127 The variations of the O-Ge-0 bond angles are larger than those of O-Si-O [Figs. 4(a) 

128 and 4(b)], whereas the tilt angles in Si0 2 vary more significantly than those in Ge0 2 

129 [Figs. 5(a) and 5(b)] with respect to the volume. These results indicate the strong 

130 rigidity of the O-Si-O bonds. The experimental study reported that the O-Ge-0 bond 

131 angles markedly change in q-Ge0 2 , while the tilt angle in q-Si0 2 varies significantly as 

132 the pressure increases, which agrees well with our result. 8 The metallic property of Ge, 

133 as mentioned above, is also attributed to the distorted Ge0 4 tetrahedra and the variation 

134 of the O-Ge-0 bond angles from the ideal tetrahedral angle. This characteristic of the O- 

135 Ge-0 bonds leads to markedly reduced lattice stress at the Ge/Ge0 2 interface during the 

136 oxidation process compared with its Si counterpart. Kageshima and Shiraishi reported 

137 that Si0 2 forms a cristobalite structure at the initial stage of Si oxidation transforms into 
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138 a quartz structure after the Si atom at the Si/Si0 2 interface has been ejected to release 

139 the interface stress. 4 We also investigated the emission of Ge atoms from the Ge/c-Ge0 2 
wo interface following their model and found that Ge atoms are hardly emitted. 5 We have 

141 concluded that this is because the dispersion of the bond angles around the suboxidized 

142 Ge atom at the Ge/Ge02 interface is larger than that around the Si atom at the Si/Si02 

143 interface. The present result that the bond angles around Ge atoms in c-Ge0 2 also more 

144 drastically change than those around Si atoms as well as those in the other phases supports 

145 the conclusion in our preceding study. 

we 4. Conclusions 

147 We have calculated the bond lengths and bond angles of the rutile, a-cristobalite, and 

148 a-quartz phases of Ge0 2 and Si0 2 using first-principles electronic-structure calculations. 

149 It was found that rutile Ge0 2 is the most stable phase among the structures of Ge0 2 

150 examined here while Si0 2 preferentially forms the quartz structure, which agrees well 

151 with previous first-principles results. Symmetry allows four different O-Ge-0 (O-Si-O) 

152 tetrahedral angles in Ge0 4 (Si0 4 ) in the case of the cristobalite and quartz phases, and 

153 the angles in Ge0 4 are more distorted from the ideal tetrahedral angle (109.5°) than those 

154 in Si0 4 at the equilibrium volume. Moreover, we have also examined the variation of the 

155 bond lengths and bond angles with respect to volume and found that the mechanisms 
we leading to compression and expansion are markedly different between Ge0 2 and Si0 2 

157 even though the volume compressibilities and expansibilities are almost identical: the 

158 tetrahedra of Ge0 4 are significantly deformed under pressure whereas the tilting angle 

159 composed of two tetrahedras markedly varies in the case of Si0 2 . These characteristics of 

160 Ge0 2 , i.e., the ground-state phase of the oxides and the difference in the variation of bond 
lei angles with respect to the volume, can be interpreted in terms of the metallic properties 

162 of the bond network of Ge. Thus, our results are highly relevant to the low defect density 

163 at the Ge/Ge0 2 interface because the deterioration of the interface is suppressed owing 

164 to the flexibility of the O-Ge-0 bond angles. 
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Tabic I. Calculated and experimental lattice constants (unit: A). 

Present work Other works Experiment 



a c a c Ref. a c Ref. 



q-Ge0 2 


4.897 


5.636 


4.870 


5.534 


20 


4.987 


5.652 


34 


c-Ge0 2 


4.818 


7.128 








4.985 


7.070 


6 


r-GeOs 


4.418 


2.886 


4.283 


2.782 


20 


4.397 


2.863 


35 


q-Si0 2 


4.850 


5.348 


4.883 


5.371 


20 


4.916 


5.405 


9 


c-Si0 2 


4.925 


6.828 


4.950 


6.909 


20 


4.929 


6.847 


36 


r-Si0 2 


4.147 


2.662 


4.175 


2.662 


20 


4.180 


2.667 


37 



Table II. Bond lengths k (in A), bond angles 6j, and tilt angles 5 (in deg) of GcC-2 and SiC-2 polymorphs. 
6i, #2, #3, and 6* 4 are the O-x-0 angles, where x represents a Ge or Si atom, k and 9j are assigned 
according to the magnitude. 

l\ I2 61 62 63 04 8 



q-Ge0 2 


Present 


work 


1.763 


1 


.755 


114.13 


110.69 


107.28 


105.39 


29.66 




Ref. 34 




1.741 


1 


.737 


113.1 


110.4 


107.7 


106.3 


26.54 


c-Ge0 2 


Present 


work 


1.760 


1 


.760 


120.69 


111.39 


109.95 


101.72 


35.64 


r-Ge0 2 


Present 


work 


1.918 


1. 


,887 


80.25 












Ref. 35 




1.903 


1 


.871 


80.2 










q-Si0 2 


Present 


work 


1.608 


1. 


.603 


110.58 


109.37 


109.23 


108.55 


17.85 




Ref. 9 




1.614 


1 


.605 


110.52 


109.24 


108.93 


108.81 


16.37 


c-Si0 2 


Present 


work 


1.604 


1 


,603 


111.46 


110.02 


109.01 


108.15 


25.41 




Ref. 10 




1.603 


1 


.603 


111.42 


109.99 


109.03 


108.20 


23.25 


r-Si0 2 


Present 


work 


1.786 


1. 


.751 


81.02 












Ref. 37 




1.810 


1 


.758 


81.35 
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235 Figure caption 

Fig. 1. Unit cells of quartz (a), cristobalite (b), and rutile (c). Black and white circles 

236 

are and Ge (Si) atoms, respectively. 
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237 Figure captions 



Fig. 2. Total energy per molecular unit (m.u.) as a function of volume for all polymorphs 
of Ge0 2 (a) and Si0 2 (b). The zero of the energy scale is rutile for Ge0 2 and quartz 
for Si0 2 . 
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239 Figure captions 

Fig. 3. Variations of x-0 bond lengths A/j in quartz (a), cristobalite (b), and rutile (c) 
structures from their equilibrium points, x represents a Ge or Si atom. All and Al 2 

240 

correspond to squares and circles, respectively. Black (white) symbols are the results 
of Ge02 (SiO"2) and lines are only eye guides. 
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241 Figure captions 

Fig. 4. Variations of 0-x-0 bond angles A8j in quartz (a), cristobalite (b), and rutile 
(c) structures from their equilibrium points, x represents a Ge or Si atom. A9i, A9 2 , 

242 A#3, and A84 correspond to squares, circles, upper triangles, and lower triangles, 
respectively. Black (white) symbols are the results of GeO"2 (SiG^) and lines are only 
eye guides. 
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243 Figure captions 

Fig. 5. Variations of tilt angles AS in quartz (a) and cristobalite (b) structures from 
their equilibrium points. AS is related to x-O-x bond angles. Variations of x-O-x 

244 bond angles Atp in rutile phases from their equilibrium points are shown in (c). x 
represents a Ge or Si atom. Black (white) symbols are the results of GeO"2 (SiO^) and 
lines are only eye guides. 
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